NPS ARCHIVE 
1966 

KAUFER, R. 



THE DYNMIC PROPERTIES OF THE FULL JOURNAL BEARING, 
CONSIDERING CAVITATION 
BY 

RICHARD A. KAUFER 

LIEUTENANT COMMANDER, UNITED STATES NAVY 



MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
MAY 1966 

_ PROFESSOR HERBERT H. RICHARDSON 

Thesi s 
K1465 






'.4 



. ' * 



I 

"•TO . 



tc School 



I 



THE DYNMIC PROPERTIES OF THE FULL JOURNAL BEARING 



CONSIDERING CAVITATION 
BY 

RICHARD A. ^AUFER 

LIEUTENANT COMMANDER , UNITED STATES NAVY 
B.S., MECHANICAL ENGINEERING, MARQUETTE UNIVERSITY (1957) 
SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS 
FOR THE DEGREE OF imAVAL ENGINEER 
AND THE DEGREE OF 

MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

AT THE 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
MAY 1966 



Signature of Author: 



Department of Naval Architecture and 
Marine Engineering, 20 May, 1966 



Certified by; 

Herbert H. Richardson 

Associate Professor of Mechanical 

Engineering 



Accepted by; 



Chairman, Departmental Committee of 
Graduate Students 



it 



,1 Scliool 

ABSTRACT 

THE DYNAMIC PROPERTIES OF THE FULL JOURNAL BEARING, 

CONSIDERING CAVITATION 
by Richard A, Kaufer 



Submitted to the Department of Naval Architecture 
and Marine Engineering and the Department of Mech- 
anical Engineering on 20 May, 1$66 In partial 
fulfillment of the requirements for the degree of 
Naval Engineer and the degree of Master of Science 
In Mechanical Engineering. 



The unsteady notion of a circumferentially grooved journal 
bearing Is analyzed and programmed on a digital computer, using 
the short bearing approximation of Reynolds* equation. The 
boundary conditions of continuity are used, and particular em- 
phasis is placed on determining the effect of the size and 
shape of the cavitation zone on the unsteady characteristics of 
the bearing. 

Tabulated data consisting of bearing loads, load angles, 
and spring and damping coefficients are presented for a wide 
range of bearing dimensions and operating conditions. A com- 
puter program Is also written which can be used to obtain a 
wider range of values than are tabulated, and/or as a sub- 
routine in an analysis of a journal bearing supported system. 

The analysis shows that the dynamic properties of a jour- 
nal bearing are greatly Influenced by the shape and extent of 
the cavitation zone, thus Indicating a need for experimental 
work with the primary object of ascertaining film behavior 
under unsteady conditions. 



Thesis Supervisor; 

Herbert H. Richardson 

Associate Professor of Mechanical Engineering 



11 



ACKNOWLEDGEMENTS 



The author expresses his appreciation to Professor Her- 
bert H, Richardson, Department of Mechanical Engineering, not 
only for his guidance and counsel but also, and in particular, 
for suggesting the subject of this thesis. 

The author also acknowledges credit to the MIT Computation 
Center through whose facilities a considerable part of this 
thesis evolved. 

Finally the author is indebted to the United States Navy 
for granting him the opportunity to pursue graduate work at 
the Massachusetts Insititute of Technology, 



ill 



TABLE OP CONTENTS 



Abstract 

Acknov/ledgements 
Table of Contents 

Nomenclature v 

Introduction 1 

Analysis of Unsteady Motion ^ 

Results 36 

Discussion of Results 

Conclusions ^9 

Appendices 50 

A. Short Bearing Approximation 51 

B. Computational Procedure $2 

C. Tabulated Data 63 

D. Table of Integrals 10? 

Bibliography 108 



iv 



NOMENCLATURE 



BNM 

c 

CB 

CJ 

CJ» 

d 

D 

e 



2F 




h 



h* 



I 

KNM 

2L 

n 



damping coefficient, as in Bj>f.; Ib-sec/in, 

B/n/?n c 

dimensionless damping coefficient, 

PdLR 

radial bearing clearance; in. 
bearing center 
journal center 
perturbed Journal center 
perturbation of Journal center; in. 

Journal diameter; in. 

Journal eccentricity; in, 
total bearing load; lb, 
total steady state bearing load; lb, 
dimensionless bearing load, 

PfiLR 

bearing forces in 1 and 2 directions respectively; lb, 

steady state bearing forces in the 1 and 2 directions 
respectively; lb, 

bearing forces in r and s directions respectively; lb, 

radial bearing clearance at any point; c(l+ncos©); 
in, 

radial bearing clearance along the upstream boundary 
of the cavltated zone; c(l + ncosO*); in. 






h*dz 



spring coefficient, as in K^.^; Ib/in. 

dimensionless spring coefficient, — 

p^LR 

total bearing length, excluding circumferential 
groove; in, 

eccentricity ratio, 

c 



V 



p 

Po 

p’ 

Pa 

PR 

R 



r,s 



u,v,w 



V 

w 



x.y 



x,y,z 



‘b* 

0 

UJ 

A* 

A 

e 

X 



pressure at any point; psla. 

feed pressure; psla. 

feed pressure, (p^^ “Pa^* pslg. 

ambient pressure; psla. 

pressure ratio, — 

Pa 

yolxmetrlc flow In the m direction; In^/sec. 

Journal radius 

radial and tangential reference axis 

fluid velocities In the x.y, and z directions res- 
pectively 

perturbation of Journal center; In/sec. 

film width at any point measured from the circumfer- 
ential groove; In. 

horizontal and vertical reference axis 

film coordinate axis 

_d^ 

dt 

w’idth of a oil strlatlon at the downstream boundary 
of the cavltated zone; In. 

width of a oil strlatlon at the upstream boundary of 

the oavitafced gone} in. 

angle between bearing load and line of centers; ra- 
dians 

rotational speed of Journal; radlans/sec. 

relative tangential velocity between Journal and 
bearing, (tu- 20); radians/sec. 

dimensionless squeeze, 

Pa 

viscosity; Ib-sec/in®. 

angle between line of centers and any point on the 
J ournal ; radians 

shear stress; psla. 
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INTRODUCTION 



Background 

The knowledge of the dynamic properties of a Journal bear- 
ing is necessary in the design stage both to ascertain the 
Inherent stability of the bearing itself and to predict the 
overall response of a bearing supported system to external and 
Internal excitations. Some of the many problems which can be 
encountered are self excited half speed whirl, excessive am- 
plitude of vibration of the Journal within the bearing near 
the critical speed of the system, and the Inability of the 
system to recover from a dynamic load. 

Since Reynolds 'derived the three dimensional, continuity, 
differential equation for determining the pressure distribu- 
tion in an oil film betxveen two moving surfaces, all Journal 
bearing investigations have used Reynolds' equation, usually 
in one of two simplified forms because of its mathematical 
formldabillty , In the two dimensional or infinite bearing sim- 
plification, the pressure Induced flow in the axial direction 
is neglected. In the short bearing simplification, the pres- 
sure Induced flovj in the circumferential direction is neglected 
by assuming that it is much less than the velocity induced flow 
caused by the relative motion between the Journal and the bear- 
ing. 

However , ‘ defining and applying the proper boundary con- 
ditions has proved to be the most difficult task in obtaining 
a good solution to Reynolds' equation, even in one of the re- 
duced forms and under steady state conditions, Sommerfeld, in 
1904, was the first to obtain an exact solution to Reynolds' 
two dimensional equation, but his solution allowed negative 
pressures in the diverging portion of the bearing. This led 
to the theory that \mder certain conditions of loading, vapor- 
ization could occur because a fluid is unable to support any 
appreciable tensile force. This theory has since been verified 
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by making observations of tranlucent bearings which clearly 
show film rupture[l] Theoretically, the cavita- 
tion pressure should be the vapor pressure of the fluid, but 
In actual practice the cavitation pressure has been found to 
be very nearly ambient pressure because of entrapped air and 
other gases. 

In 1952 , 0ckvlrk[5] developed a steady state solution of 
Reynolds’ equation using the short bearing approximation and 
asstimlng that the entire unloaded half of the bearing vjas ca- 
vltated and at zero pressure. More recently, Pedor[63 and 
Donaldson[7] obtained solutions of Reynolds' three dimensional 
equation using series solutions and applying Sommerfeld's 
boundary conditions, but requiring that the feed pressure be 
sufficiently high so as to suppress cavitation. However, this 
Is an imreallstic situation In most bearing applications, par- 
ticularly those in vrhlch stability is a problem. 

Two recent analyses of unsteady motion are of special 
Interest in this investigation. Holmes^S] extended Oclcvirk's 
solution to the particular case of a rigid shaft supported by 
S 3 nimietrlcal bearings, Rlchardson[9] also used the short bear- 
ing approximation, but defined the cavitation boundary more 
rigorously by assuming that there was a zero pressure gradient 
along the boundary and by assuming that the shape and extent of 
the cavity was a function of feed pressure and squeeze , -Rich- 
ardson obtains a solution to his analysis by digital methods. 
However, both Holmes and Richardson violate continuity since 
they do not consider the oil entering and subsequently leaving 
the cavitation zone. 

In 1957# theories which took Into account continuity of 
flow through the cavitation zone were Independently developed 
by Floberg[l03 and Rightmlrej^ll] . Experimental work by Flo- 
berg and Ishll[4] has Indicated a close agreement with the 
continuity theory under steady state conditions. 
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Purpose 

The purpose of this Investigation Is to analyze the un- 
steady motion of a circular, circumferentially grooved journal 
bearing, using the short bearing approximation and the bound- 
ary conditions of continuity. A circumferentially grooved 
Journal bearing was chosen because its line source makes the 
feeding arrangement independent of the journal attitude within 
the bearing. The results of this analysis will be presented 
In two forms: a computer program which can be used as a sub- 
routine in a systems analysis, and In the form of tabulated 
data. of spring and damping coeffeclents for use In hand calcu- 
lations , 
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ANALYSIS 0? UNSTEADY MOTION 



As s urn •nt ions 

In arriving at a solution of Reynolds* equation which 
v;ill adequately describe the unsteady motion and dynamic pro- 
perties of a circumferentially grooved Journal bearing, the 
follov/ing assumptions will be made: 

1. The fluid is Newtonian and Incompressible. 

2. The viscosity is constant. 

3. The flow is laminar. 

L. The mass of the fluid is negligible. 

5. Surface tension' is negligible. 

6. The fluid adheres to the Journal and bearing surfaces, 

7. There is no distortion of the Journal or bearing 
surfaces due to hydrodynamic pressure; only light 
and moderately loaded bearings will be considered. 

8. The fluid film is so thin that it has no curvature. 

9. The fluid film is so thin that there is no pressure 
gradient across the film. 

10. The vapor pressure of the fluid is equal to ambient 
pressure . 

11. The pressure induced flow in the circumferential di- 
rection is negligible compared to the velocity 
induced flow in the circumferential direction 

12. In the cavitated region, the fluid divides into an 
Infinite number of streamers. 

The last two assumptions will be discussed in greater 
detail in later sections; the remaining assumptions are self 
explanatory . 

Journal Bearing Oil Film Theory 

Figure I illustrates the general case of the unsteady 
motion of a Journal within a fixed rigid bearing. At any 
instant in time the Journal has a rotational speed CU about 
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Figure I 

Joxurnal -bearing motion referred, to a fixed point in space 




Figure II 



Journal -bearing motion referred to the line of centers 
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Figure III 

Oil flla forces in the x direction 




Figure IV 

Oil film forces in the z direction 



P + 



iL 




% 

its own center C J , an orbital speed of 0 about the bearing 
center CB and a radial velocity e. However in analyzing this 
systera, it is raore convenient to fix the line of centers in 

space and use that line as a reference as shovrn in Figure II. 

» 

Then the journal has a clockwise rotational speed of (UJ - 0) 
and the bearing has counterclockwise rotational speed of 0. 

The forces acting in the x and z directions on an ele- 
mental cube dx dy dz are shown in Figures III and IV, where 
u, V, and w are the fluid velocities in the x, y, and z di- 
rections respectively au any point. Equilibrium of forces 
in the X and y direction yields; 

dx dy dz 



dP _ d'Cz . 
dz by ^ dx 



V/ith the assumption that the film is so thin that has no 
pressure gradient across it; 



^ - 
dy 



0 



The sheer stresses are obtained from Newton's law for laminar 
viscous flow[l2]: 






< ^ 
dy 




dw 

dy 



Substituting these expressions into equations 1 and 2 gives: 
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( 3 ) 



^ +yj .^liL _ feu , ^^;ul 



rtP ..^.“^w 



=// 



.x^)z 



6 

^y 



N \ 

evx , gw 

<V 1 



(4) 



Asain considering the i'act that the film is very thin in the 
y direction, so that the velocity gradients in the y direc- 
tion are much greater than the velocity gradients in the x 
and z directions: 

\ 

d W dVX 

by 

bu bu 

bz jjy 



Then equations 3 and 4 become: 



bp _ w b^u 


(5) 


p P = yt/ p ^ W 

bz ^ W 


(6) 



With the assumption of constant viscosity, equations 5 and 6 
can be integrated twice: 



u 



1 ^ 
27? bx 



+ C,y + C. 



(7) 



w 



-i- v" 

2/f ;^z ^ 



+ Cjy + 



( 8 ) 
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where C, , C-^ , and C*, are constants of Integration. 

From Figure II, the boundary conditions are: 

aty = 0: u = (CO - 0)R, v;= 0 

at y = h : u = - , w = 0 

where h is the film thickness at any x. 

Applying these boundary conditions to equations 7 and 8: 



With the assumption that the pressure induced flow in the 
circumferential direction is much less than the velocity in- 
duced flow due to the relative motion between the journal 
and the bearing, equation 9 becomes: 



The boundary Velocities in the y direction are: 

at y = 0 , V = y 
at y = h, V = 0 

By continuity: 

ou ^ ^ ^ 

by bz 

Substituting equations 10 and 11 into the continuity equation 




( 9 ) 




( 10 ) 



u = -^y + (CA)- 0)R 



( 11 ) 



9 



A 



0 



h 



y+ (LL) - 0)Ridyi + ~{y dy) 



J 



*’y>j 1 = 



Integrating the above expression from y = 0 to y = h and 
applying the boundary velocities in the y direction, the 
final result is the short bearing approximation of Reynolds* 
equation for unsteady motion: 



^z3 ^x " h" ^ 



(12) 



where -A_ = (LU - 20) 



This equation is much easier to manipulate and solve 
than Reynolds’ complete equation, and for moderately loaded 
bearings with length to diameter ratios of less than | the 
loss of accuracy introduced by the short bearing simplifi- 
cation can be expected to be of the same order of magnitude 
as some of the other assumptions which v:ere necessary to 
make. A comparison between an exact solution obtained by 
Donaldson ([73 under steady state conditions (0=0, e = 0) 
and the results of this investigation is presented in 
Appendix A. It is proposed that the same deviations could 
be expected in the case of unsteady motion between exact 
and approximate solutions. 

Again assuming constant viscosity, equation 12 can be 
Integrated twice with respect to z to obtain the pressure 
distribution on the bearing surface: 



3/;R-A. Vo SJ-f . 



+ Cj z "f C 2 



(13) 
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where C, and are constants of Integration to be de- 
termined by the film boundary conditions. 

Expressions for the unit flovrs in the x and z directions 
will be required in the next section for determining the cav- 
tation boundaries. These expressions can readily be obtained 
by integrating the velocities u and v across the film as 
follows : 

Q^. = I u dy 

“'O 

r" 

Q- =i V dy 
-*0 



Substituting in equations 10 and 11 the final result is: 



Qx 



_hf_ ^ 
12/1 bx 



Q:, = (03- 20) ^ 



h“ ^ 
l^/■i bz 



(14) 

(15) 



Film Boundary Conditions 

The unwrapped oil film for one half of a circumferen- 
tially grooved journal bearing is shown in Figure V. The 
approximate shape of the cavitation zone as observed experi- 
mentally [l], [ 2 ], [3l» [ 43 , is defined by the boundary 
1-2-3; steady motion, point 1 occurs immediately 

downstream from the point of minimum film thickness (x = ttR) 
where the rate of film expansion is the greatest, but in the 
case of unsteady motion it must be assumed that cavitation 
could commence at any x. Since the vapor pressure of the 
fluid is assumed to be ambient, then the entire cavity must 
be at ambient pressure, with the result that one of the 
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.1 





boundary conditions is: 



at z = vr, p = Pa 

where v; Is the film vjldth at any x and is equal to L 
where cavitation is not present. 



Applying this boundary condition to equation 13, the result is 

(16) 



[3 6//.1 2 ax ^ / 

P - Pa=| ^ ) + C, (z - w) 



^x h-^ 



Fiprure V 

The unwrapped oil film 

^ 1 3 4 




z 

ttD >- 

>- x=R© 

With the assximption that the fluid adheres to the sur- 
faces of the journal and bearing, the fluid must split into 
one or more striations upon reaching the boundary 1 - 2 in 
order to satisfy contlnuitjr and to prevent subambient pres- 
sures from occurring. This is shown in Figure VI for the 
particular case where h^>h, . 
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Figure VI 

Continuity of flow across the cavitation boundary 




The upstream boundary 1 - 2 is determined by considering 
the flov7 through an elemental volume h dx dz at A. Then the 
flow entering the boundary at A Is : 

From equations l4 and 15 this becomes: 

Q,n = 20,)Rhdx-^|^ ^ dz 

Since there are no pressure gradients in the cavitated zone, 
the only flow is velocity induced flow caused by the rela- 
tive motion betxveen journal and bearing. Thus the flow 
leaving the boundary at A must be : 

= 1(0) - 2^)Rh dx 

Therefore, by continuity of flow across the boundary: 



13 



p 

z 



0 



A 



or : 



P = 0 



z 



w 

along boundary 1-2. 

Applying this boundary condition to equation l6 to evaluate 
the constant C, the pressure distribution in the region 
bounded by z = 0 and 3-'^“l-2is found to be: 



where p * = p^ - p 

The first term on the right side of equation 1? is the 
hydrodynamic pressure distribution and the second term is the 
hydrostatic pressure distribution. The slopes of these two 
terms at z = w, are, respectively: 



The existence of cavitation at any x upstream from point 2 and 
downstream from point 3 can be ascertained by letting w = L in 
the above expressions. Then the condition for cavitation to 
occur is : 




(17) 




( 18 ) 




(19) 




(20) 
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i.e. the positive slope of the hydrodynamic pressure Is 
greater than the negative slope of the hydrostatic pressure. 
This is Illustrated In the left diagram of Figure VII. 

Figure VII 

Axial pressure distribution In region 1-2 





If the condition for cavitation is satisfled(equation 20) then 
w adjusts so that: 



\ 


<»p' 


/hd 





0 



Substituting equations 18 and 19 into the above expression, 
the film x^ldth along boundary 1 - 2 is : 



= 



p'h^ 



3>^(R^|~-2y) 



( 21 ) 
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Making the substitutions: 




y 



e cos © 



Equation 21 becomes: 




3//(cn^sin©~ 2e cos 6) 



( 22 ) 



If equation 20 is not satisfied, then: 
w = L 

It can be seen that equation 22 is invalid for the down- 
stream boundary 2 - 3 t>y considering the steady state case. 
Then equation 22 becomes: 



Under these conditions, the expression under the square 
root becomes negative at angles greater than 0 = 2rr, implying 
that the cavitation zone cannot extend past 6 = 2n, Physically 
this would not necessarily be true for high journal speeds and 
low feed pressures. 

The downstream boundary is determined by considering con- 
tinuity of flow through the cavity as illustrated in Figure . 
VIII. Let the width of one striation at the upstream bound- 
ary be Az, Then the amount of fluid in the striation as was 
previously determined is: 

IvV.Rh'*^ A z 




cncu sin 6 
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where h* Is the film thickness at z = w on the upstream 
boundary . 



Figure VIII 

Continuity of flow across the dovmstream boundary 




V/ith the assumption that the fluid adheres to the Journ- 
al and bearing surfaces, then by continuity the same amount of 
fluid must leave the cavity at z = w at point B on the down- 
stream boundary. Since we have assumed that there are no 
pre§surg gradign'^s in the gavitjr, '^here can bg no raoljion Qf 
the strlation in the z direction. 

Summing up the total flow in the unit volume h Ax Az 
at B, the result is : 



|ARh^:*A z 



ill hR 

12M hz 



Ax 



V/ 



From eq\jatlon 17: 



|ARh Az 



(24) 



bz 



iw 



3M 

h”^ 



(rA 



bx 



2y)w - 



£i 

w 



i 



(25) 
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Substituting equation 25 into 24: 



lARh^Az - ^(rA— - 2y)w - ^ Ax = lARhAz 

12/f ^x w 



Ass\imlng an infinite nximber of strlations, and then taking the 
limit as Ax— ^dx and Az— <»'dz: 



|ARh* dz - 



|(ha|| - 2y) 



h^P* 

12/y w_ 



dx = 



|ARh dz 



( 26 ) 



Making the substitutions 



dx = R de 



dz 

6x 



= dw 



w 



on 

R 



sin 6 



y = - e cos © 



Equation 26 becomes 



•|ARhMw - 



^w(- cnAsin 6 + 2e cos 6) 



_ h^p« ] 

12 wj' 



R de 



= l-y^-Rh dw 



or; 

dw _ 5 w(cnAsin 6 - 2e cos 6) + 6// vr 

de ~ , f. 371 

A. ( h “ h'*'" ) 



(27) 



The film width, w, in equation 27 is in terms of two vari- 
ables, dw/de and h*, and so must be solved by an iteration 
method. 
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Oil Film Forces 



The oil film forces which vJill be considred are: 

1. Hydrodynaraic , hydrostatic and shear in the complete 
film region bounded by z = 0 and 4-l-2-3-'^ 
in Figure V. 

2. Shear in the cavitation zone. 

The differential hydrodynamic and hydrostatic force on a 
strip of the journal w, RdOis: 



d-Fe. = / p dz R d0 

i 

Jz-c 

Substituting in p from equation 1?: 

Cr'' 

= I I - 2 y)(z^ - z»,0 






p'd - -) I dzRd© 



which when integrated becomes: 



=! - 



w 






''h 






2h 



•(RA|^ - 2y) -i- |-p'w 



R de 



Resolving this differential force in the 1 and 2 directions as 
shown in Figure IX: 



dFp, = dFy cos 6 
r< r 

dFp^ = dpj. sin e 

The tangential shear stress in the complete film region 
due to the relative sliding velocity is (see Figure II): 

Xi = "^(CUR - ‘0e cos 0 + e sin ©) 
h 

However, ^e^CUR and e«(OR, so that: 

•T- ^ A^'VR 
~ h 
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Figure IX 

Coordinate systen for oil film forces 




Then the differential force due to shear on a strip w, Rd©is 



dF- = 'C wR d© 

w 



Substituting in the value forX»: 



■o h 



Resolving the shear force in the 1 and 2 directions: 
dF„, = - dR-sinGde 
dPr^ = dF,- cos & d© 

In the cavitated zone, there is no hydrostatic pressure, 
since p = p^ . Also, since there are no pressure gradients, 
there is no pressure flow and thus no hydrodynamic pressure 
exists . 
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The shear forces in the cavltated zone are obtained by 
assuming an infinite number of striations as in the preceding 
section. Refer to Figure X v;here one striation is shovm. 



Figure X 

Continuity of a single oil striation 




By continuity: 

|ARh*S^ = iARhS 

where and S are the striation widths at the upstream 
and downstream boundaries respectively. 

Then: 

h*S-"- = hS 

or : 




The differential shear force due to one striation at any 6 is 
dfr ='CSRde 

o 
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Substituting in the expression for 8 fron above: 



df. ^Rd0 



h 



h 



In the limit, as S^''—>dz: 
df<- = — — h-dz d© 



Then the differential force due to shear across the complete 
vridth of the cavitated zone at any 0 is: 

I h® 

"'z » VV 



h* dz d© 



vrhere the integral must be evaluated numerically. 



Resolving the cavitated zone shear force in the 1 and 2 direc- 
tions : 



dFc =-dFr sin©d0 

“Cl 

dFc = dFr cos 0 d© 
Sc2 •'C 



Susmins up the- forggs in the i 2 dlrgstlsns gM m&lclng th§ 

substitutions; 



6^ 

'dx 



- ^sine 
R 



y = - e cos 0 



Then the differential forces in the 1 and 2 directions are: 



dF, = 



w^R 



2h" 



(,/Lcn sin 0 - 2e cos 0)cos 6 + — — cos 0 






h 



h‘ 



d© (28) 
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sin 6 



dP„ = 

z 






2 h" 



(-A,cnsine - 2ecose)sin6+ ~ ■■ 



x^rhere I = j h*dz 
Jw 



Then: 



/yiVn^vT AiiUR^I 1 

+ — ; cos © + r — COS G d© 

h J 



r2,r 
F, =/ dP, 

P =/ dF^ 

Jo 

which must be evaluated numerically. 



( 29 ) 



The attitude angle is 
0 = tan ' (^) 



( 30 ) 



Spring and Damping Coeffecients 

Refer to Figure XI. At a given instant of time the jour- 
nal center is at CJ and is being acted upon by the forces F, 
and p2 which were evaluated in the previous section. Now the 
journal is given a perturbed displacement, d, and a perturbed 
velocity, V, to the position CJ’; then the forces acting upon 
the journal are of the form: 

Fj. = F,p +Kj.(r,s)d+ Bj.(r,s) 

Pf^ ~ ^2,0 "^^5(^»^) d + B^(r,s) 

where F,j> and P^^are the steady state forces at C J ; and K 
and B are spring and damping coefficients. 
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Figure XI 

Perturbation of journal center 



C 3 




In evaluating the previous equations it is first nec- 
essary to transfer the forces ?, and to the r and s axis. 
However, it is more convenient to continue using differential 
forces; thus: 

dFj. = dF, cosc<._ dF^sinc^ 

dFg = dF.^sin»-A+ dF^cosc^ 

Since a small perturbation is being considered: 

CJ CJ» « CB CJ* 



Then: 



cos 0< 1 
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with these approrcinatlons , the differential radial and tan- 
gential forces become : 



dF^ = 


df, - dPj 5k 


(31) 


dF, = 


1?, ^ ^ -iPl 


(32) 



Substituting the eo^uations 28 and 29 into equations 31 and 32: 

r 



dF^, = ^ — cos © 

I 2 h-^ 



- 2e cos 6 + yV nc sin 6 



+ I p’wR cos 0 

3lne sine - ^ Cosine 

h h** cn 2h^ 

r . ^ 1 

- 2e cos 0 + _/'i- nc sin © | + -g- p’wR sin 6 



I 



, /^0)R®w ^ I 

-f cos 6 -j- r — COS 0 / > d© 

h J ) 

dFr. = ^ cos 0 I - 2e cos © + -/Vcn sin 6 

s I cn \^2h^ [ 



+ g’p'wr cos © 



//a.’R 2 w , ^ /fw^R , ^ 

• — ^ii-Jlsin© - — Sin ©j + Sin© 

h h^ / 2 h^ 

— 2e cos 0 + vL nc sin ©1 + ^ p’wR sin © 



MiX)R^ 



w 



h 



cos © + 



/yaiR^i 



h^ 



cos ©^ d© 



Let : 



• « 

e = r 

h = c (1 + n cos ©) 




Substituting the above into the equations for dF^ and dF^ 

I •«» 
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and rearranging: 



dF^ = 



, >/ R r 

cos e 



2r cos 6 



iO nc sin 0 



2q 3 ^ (1-t-ncosG)^ (1 + ncosG)'^ 

2s sin 6 



1 3 , - ^ R^w sin 6 

r rr- !+ tp'wRcOS© - 7- -^7- 

(1 + n cos 6)“j (lT-ncos6)c 

LU sin 6 J-l vr^R sin 6 
c^ (1 -r- n cos 0)^ 

(,-L's sin 6 



-2c^ 

2ss sin 0 



2rs 



(l + ncos0)^ end -f n cos 6)^ j 

R^ws cos 6 a/LUR^Is cos 6 ) 

. \ do 



end + n cos 6)^ 
1 p?wRs sin 0 



2cn 



dF. = 



c^n(l + n cos ©) c 


^n(l + n cos 6f j 


o 

o 

CO 

CD 


2rs cos © 


2c^ 1, (l+ncos©)^cn 


2ss sin 0 


p’ swR cos 6 

• _ . 


end + n cos 0)^ 


2cn 


MW R^ls sin 0 


M w^R sin © 1 

HIM .J- - ■ 1 


c^n(l + n cos 0)^ 


2c^ i 


iU nc sin 6 


2s sin © *j 


(1 + n cos 0)^ 


(1 -f n cos 0)^ J 



(33) 



LU s sin 0 
(1 + n cos 0)^ 

MUJR^vrs sin Q 
c^nd + n cos 0) 

2r cos © 

(1 + n cos 0)^ 

+ ^ p?wR sin © 



cos 6 /-flUR^lcos e 
+ -{- 



c(l+ncos0 c^(l+ncos0)^ 



-d© 



(34) 



Equations 33 and 34 are non linear in n. However, with the 
assxmption that the perturbation is small, the equations can 
be linearized. From Figure XI: 



cn^ + r = cn cos of. cs on 



or : 



n 



n<> + 



r 

c 
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where — « n^and Is time dependent. 

By Taylor's theorem: 

= f(n.)+ |f'(n„) 

X* 

neglecting second order and higher poxvers of — . 

c 



Then, by also neglecting other second order terms: 



r f (n^ -^l-) ^ 

I +|) I f(nc ) 

s f (n^. + |-) ^ s f (Uq ) 



r s = 



0 



s s = 0 



Linearizing the terms of equations 33 and 3^ as above: 



dP^ = 



// vr^R 
2c^ 



cos e 



2r cos © 



UJ cn sin © 



(1+ncos©)^ (1+ncos©)^ 



4 - sin 0(1+ n cos © - 3noos©) 2s sin Q 

( 1 + n cos (l+ncos©)'^ 



+ p’wR cos © - 

MUJR^I sin e 
c® (1 + n cos ©)^ 



y^HOR^vi sin © 



AUVR^viv sin © cos © 



c(l+ncos©) c^(l + ncos©) 

2MUJR^lra sin 6 cos 6 

c^ (1 + n cos ©)^ 
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dp., = 



■^1 w^Rs66>sin^ 6 
2c^ (1 + n cos 6)^ 

MLU R^Is cos 6 
c^n 0- + n cos 6)^ 



p^wRs sin © // 6(/' R^ws cos 6 



2cn 



• d© 



c^n 0- + n cos 6) 

(35) 



y^'/ w^RstUcos 6 sin ©_^ p^svjR cos 6 MtU^ v:s sin 0 



2c"^ (1 + n cos 

yq U) R^Is sin © 
c^n (1 T n cos 6)‘ 



2cn 



ylf w^R sin © 

2^^ 



c^n ( 1 + n cos 0) 



2r cos © 



(1 + n cos ©)' 



CU n c sin © sin © ( 1 + n cos © - 3n cos © ) 

+ — 4 - 



(1 T n cos ©)^ 

2s sin © "l 
(1 + n cos ©)^ j 

yWW/R^wr cos^ 6 



(1 + n cos ©)'* 

//40 ■/!' cos © 



+ § p ' wR sin 0 + 



cos © 



c ( 1 + n cos ©) 



c^(l + ncos©)^ c^(l+ncos0)‘ 



RA'O'R^Ir cos^ o') 



c^ (l + n cos ©)^ 



1 



d© 



(36) 



Equations 35 and 36 are of the form; 



p> 

II 


= - dKj.^r 


- dB^^r 


- dK,,s 


+ dB^jjS 




dF., 5 


: = diC . § 
s*s "" 


= dB,, s 
— - 


+ dK.. . r 
- 


1= 





where : 



1Krr = 



>/ Lb w^ R 



2c^ (1 + n cos ©) 
MUOR^v^ sin © cos 6 



sin 0 cos 0( 1 + n cos © - 3 j^cos 0) 
2//i^^’R^I sin 0 cos 6 



d 3 ^ = 
tr 



c^ (1 + n cos 0)^ 
w^R cos^ 0 



c^ (1 + n cos 0)' 



■d© 



(1 n cos ©)' 



d0 



(37) 

(38) 
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dK 



ts 



^^5 



dK. 



sin^ e 

' 2 c^ (1 + n cos ©)' 

M Uy R^I cos 6 
c^n (1 + n cos &)^J 

J-/y:^R cos 0 sin 6 
c*^ (1 + n cos G)*^ 



p’wRsine J/UJR^w oos e 

+ — + 



2 cn 



c^n ( 1 + n cos ©) 






d© 



d© 



w^RiO cos © sin © p?wR cos © /^LU R^v; sin 6 

+ 



2 c^ (1 + n cos ©)^ 
yi-flV R^I sin e 



2 cn 



c®n (1 + n cos 









c^n (1 + n cos 6 )^ 
^ w^R sin^ 6 



d© 



c^ (1 + n cos 6) 

>/g> X'j^R sin^ © 
2c^ (1 + n cos 6) 

/■/iV R^vx cos^ © 



^ d© 



(1 + n cos © - 3 n cos ©) 
ZM'JUR.^1 cos^ © 



d 3 



- 



c^ (l + n cos 0 )^ c^(l+ncos©)' 

J'i w^R sin © cos © 

. d© 

c"^ (1 + n cos ©)^ 



d© 



then: 




which must be evaluated numerically. 



The final result is 



Pv^=-K r-B r-K^s+B s+F 
n »>x- ts rs 



(39) 

( 40 ) 

e") 

( 41 ) 

( 42 ) 

(43) 



(45) 



29 



( 46 ) 



Fc = - s - B..S+K. r+B r + 
s ss iS i>t^ sr 2(7 

However, it is generally more convenient to work with 
cartesian rather than polar coordinates. Therefore the in- 
tegrated forces, and spring and damping coefficients will be 
transferred to the x - y axis as shown in Figure XII. 

Figure XII 

Journal motion in x - y coordinates 




r = y sin 0 +xcos0 


( 4 ?) 


s = y cos 0 - X sin 0 


( 48 ) 


• • • 

r = y sin 0 + x cos 0 


( 49 ) 


s = y cos 0 - X sin 0 


( 50 ) 


Fj^ = F^ cos 0 - Fg sin 0 


( 51 ) 


Fy= FySln 0 + F^ cos 0 


( 52 ) 
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Equations 4? through 52 are substituted Into equations 45 and 
46; 

Pj, = - (y sin 0 + X cos 0 ) - B j.^(y sin 0 + x cos 0 )+ F,^ 

- Ky.^ (y cos 0 - X sin 0 ) + (y cos 0 - x sin 0 ) (53) 

Fj. e - (y ooa 0 - X sin 0 ) - B (y cos 0 - x sin 0 ) + 

+ Kgj, (y sin 0 + X cos 0 ) + B^y. (y sin 0 + x cos 0 ) (5^) 

Substituting equations 53 and 54 into equations 51 and 52: 

sin 0 cos 0 y - Ky,.^cos^ 0 x - Bj-y. sin 0 cos 0 y 

- By-j- cos^ 0 X - Ky >3 cos^ 0 y + Ky -3 sin 0 cos 0 x 

0 • 

+ B y -3 cos^ 0 y - Bp 3 sin 0 cos 0 x + P,^ cos 0 
+ Kgjj cos 0 sin 0 y - K sin® 0 x + B 5 ^ cos 0 sin 0 y 

- B sin® 0 X - Kg), sin® 0 y - sin 0 cos 0 x 

- B^,, sin® 0 y - Bj.^ sin 0 cos 0 x - sin 0 

Fy = - Kyj. sin® 0 y - sin 0 cos 0 x - B^.^ sin® 0 y 

- By^ sin 0 cos 0 X - sin 0 cos 0 y + sin® 0 x 

+ Bj. sin 0 cos 0 y - B, . sin® 0 x + F, sin 0 

0 

- K 33 COS® 0 y + Kjj sin0 cos 0 x - Bj^ cos® 0 y 

+ B 33 sin 0 cos 0 X + Kjp sin 0 cos 0 y + K^j, cos® 0 x 

• * 

+ Bjj, sin 0 cos 0 y + Bj^. cos® 0 x + cos 0 
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Then: 



Fx “ “ y - Fo 

Fy = - Ky^y - Byyy + X - x 
VJhere : 



(55) 

(56) 



K^^cos=»0 + 

Kyy = Ksj COS^ 0 + 

Kxy= K^s cos^ 0 + 
Kyx “ ^sr cos^ 0 + 
Bxx = 0 + 

Byy= B 22 cos^ 0 - 
®xy = Bsr Sin® 0 + 



(K,.,.- Ky.^.) sin0 
(Ky.g - Kjy.) sin0 
(K^y.- Kjs ) sin0 
(Kgs - Kyy ) sin 0 
(B^j+ Bjp ) sin0 
(Byj. + Bjy. ) sin 0 
(By^y** ) cos 0 



cos 0 + K ^5 sin® 0 
cos 0 + Ky,y, sin® 0 
cos 0 + Ksy> sin® 0 
cos 0 + Kyg sin® 0 
cos 0 + Bs 5 sin® 0 
cos 0 + By^ sin® 0 
sin 0 - Bp -5 cos® 0 



(57) 

(58) 

(59) 

( 60 ) 
( 61 ) 
( 62 ) 

(63) 

(64) 



Computational Procedure 

The general case of \msteady motion of a circular cir- 
cumferentially grooved Journal bearing has now been analyzed 
and the differential forces and differential spring and damp- 
ing coefficients are in forms which can be solved numerically. 
A detailed procedure, using a high speed computer, is given in 
Appendix B, Hbvrever, the general computational procedure is as 
follows. For a given bearing geometry (L, R, C, and/y) and 
given operating conditions (CU , n and e); 
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D 



1. Determine direction of flow (± 0) and locate up- 
stream and downstream cavitation boundaries. 

2. For each 6, calculate the film width. 

3. For each 0, determine the differential forces and 
differential spring and damping constants. 

4. Integrate numerically to obtain the radial and tan- 
gential forces and spring and damping coefficients. 

5. Calculate the bearing load, load angl^ and horizontal 
and vertical spring and damping coefficients. 



Simplified Analysis 

It should be noted that equations 37 through 44, and, 
consequently equations 57 through 64 are independent of e,0, 
^e, and due to the fact that (see Figure XI): 



dK 



mm 



= a4"o k 

C7 



dB 



mm 



= A (dB^ v) 

V— ^0 ^V mm 

CX 



Since dK^„^^ and dB^^„ are functions of x-/, then w must also be 
independent of 0 and e in the above mentioned equations. 

If now frictional effects are neglected on the assump- 
tion that they are small compared to squeeze effects, then -vi- 
and CU can be used interchangeably throughout the entire anal- 
sis by virtue of the choice of reference axis. Then the equa- 
tions for the spring and damping coefficients and the equations 
for the steady state forces F^^ and F^^^ can be non-dimension- 
alized and expressed in terms of just three independent vari- 
ables; eccentricity ratio, n, pressure ratio, PR, and dimen- 
sionless squeeze , _/\f. Let: 



n = 



e 

c 
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PR = 



p 



= 



Ln 3 






Then the pertinent dlnenslonless equations are: 



/- (PR - 1) (Iv noose)' 

V 






sin 6 



Vv/ dw / VT \ a 



1 



db 




dP 

10 


, X 


p LR “ 


2 


d? 2 > 


X 


p LR “ 


• 


dK;.,C 




p LR 




dj(C pr G 


r 


p LR 


=j 


dK^j c 


1 


Pi 


~i 


dKjrC 




p LR 




d3 y* y* G 




p LR 




d3^S c 




p LR 



(1 4 - n QOS (PR - 1) 
6 nyu'"'(cos © - cos 6 '") 



( 65 ) 



( 66 ) 



L (1 -f- n cos 0 )^ L 

.T/T. , sln^ © ,w 



- -s- (i-) (PR - 1) sin© d© ( 68 ) 



‘L"^ (1 4 - n cos e)^ ■ L 



-A5 sin© cos 6 ^ ^ 

-- — (1 -}- ncos© - 351 cos ©)d© 

2 L' (1 n cos 0) 



-^^^(PR - 1 ) (g) sinejd© ( 70 ) 



■ w,, sin 6 cos 6 i , ,, ,vr. ], 

(g)^ 7 — - 7 - ^ “(PH - 1 ) (^) cosed© 

2 L (1 + ncos ©) 2 n l i 

( 71 ) 



(Wj3 sin- e 



2 P- (1 -i- n cos 6 ) 

= A 



(1 + ncos© - 351 cos©) ( 72 ) 



cos^ 6 



= A'"' 



(1 -i- n cos ©)- 
sin^ 6 

(1 + n cos ©)- 



d© 



d© 



( 73 ) 



{ 7 ^) 
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sin © cos 0 



d© 



( 75 ) 



dS r p c 
P LR 



= 



(1 V n cos 0)*^ 



dB^. = d3,^ 



The general computational procedure is the same as before, 
except that all frictional effects are neglected. 
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RESULTS 



Figures XIII, XIV, XV, XVI and XVII illustrate the shape 
and extent of the cavitation zone for various coQblnations of 
Journal motion obtained by this investigation. Results obtained 
by Ishli and Richardson ( the latter in dashed lines) are also 
presented for comparison. The direction of the circumferential 
flo/x in Figure XVI is the reverse of the circumferential flex'/ 
in Figure XV, thus Flgrre XVI is a mirror image of Figure XV. 

Tables 1 and 2 sho’.'/ a breakdox'/n of the effects of circum- 
ferential squeeze, radial squeeze, feed pressure, complete film 
friction, and cavitation zone friction on the forces and spring 
and damping coefficients for the examples in Figures XV and XIV. 

Tabulated data for the simplified analysis x/hich neglected 
frictional effets is given in Appendix C for numerous combina- 
tions of dimensionless squeeze, pressure ratio, and eccentrici- 
ty ratio. The parameters xiere chosen so as to cover as -wide a 
range of bearing applications as possible. 

Table 3 is a comparison betv/een the data in Appendix C 

r 1 

and data obtained by Holmes 8 1 who analyzed the vibration of 

L w/ 

a rigid shaft by extending Ockvirk^s analysis and boundary con- 
ditions, The data in Appendix C xcas converted to Holmes’ nomen- 
clature as follows; 



c 

2F’ 






KNM- 






2P 



2p, ""nm 



p^LR 

2F« 



V/here : 



F’ 



r 

* ( 1 - n^ ) l6n^ 

4(1 -n^)® 



1 



( 

J 



1 

2 






6 



vrhich is a fora of Ockvirk’s load nunber. The factor of 
2 is necessary because this investigation uses a circum- 
ferentially grooved bearing with a total length of 2L. 

Since Ockvirk’s theory defines the attitude angle for a given 
eccentricity ratio, it was necessary to fix both 0 and n in 
the comparison. 

The program in Appendix. B is v-;ritten in a manner so that 
any or all of the values can be called without any major changes. 
Frictional effects can be neglected simply by deleting the in- 
dicated cards, thus reducing the program to three dimensionless 
variables as in Appendix C. It should be noted that for sim- 
plic ty only one half of the bearing was considered in the 
analysis, but that the program and tabulated ^..ata are in terms 
of the entire bearing length 2L. 
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Figure XIII 
Cavitation Boundaries 
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Figure XIV 
Cavitation Boundrvje 
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Figure XV 

Cavitation Boundaries 
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Figure XVI 

Cavitation Boundaries 
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Figure XVII 
Cavitation Boundaries 
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Table 1; Effect of squeeze, pressure au:l friction on the forces and spring and 
damping constants of Figure XIV. 
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Table 2: Effect of squeeze, pressure end friction on the forces and spring and 
damping constants of Figure XV. 
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Table 3: Comparison of Koines' results[8] and tb.e results of this investigation. 
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DISCUSSION OP RESULTS 



For the steady state example of Figure XIII, Ishll ’ s results 
are very close to the results of this investigation, consider- 
ing that vjhile Ishii used continuity of flow through the cavity, 
he averaged the flovx leaving the boundary 1-2 and then used 
this average to determine the boundary 2-3. Thus it should 
be expected that Ishii ’s point C should be outside the bound- 
ary of this investigation since the average flow at vr is less 

c 

than the straight line flow, while point D should be inside the 
boundary of this investigation since the avearage flow at 
is greater than the straight line flov?. However, h does not 
vary appreciably betx'/een point 1 and 2 due to the large slope 
of the upstream boundary, and thus the average flow through- 
out the cavity is very close to the exact flow. Also, some 
minor disagreement could be expected betvjeen machine and hand 
calculations , 

The deviation between this theory and Richardson’s is 
quite pronounced along the dovmstream boundary, particularly 
in the steady state cases of Figure XIII and XIV; closer agree- 
ment is reached as the ratio of circumferential flow to axial 
flovx is reduced as in Figure XV and XVI until there is perfect 
agreement in Figure XVII where only axial flow due to radial 
velocity e is present. None of the sample runs indicated any 
differences along the upstream boundary; this vxas to be ex- 
pected since both theories use the same method for calculating 
the upstream boundary. The values for load and attitude angle 
are surprisingly close in Figure XVII, considering that dif- 
ferent numerical Integration techniques were employed for the 
tvxo theories. Thus any deviations in the results can be attri- 
buted directly’ to the differences between the doxmstream cavi- 
tation boundaries. This is best illustrated in the steady state 
examples of Figures XIII and XIV, where the most pronounced 
difference is in the fourth quadrant i ' (.1 . 5 tt <©< 2rr ) .' , 
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Richardson’s theory yields a greater film extent in this area 
and consequently a greater film force in the direction of the 
applied load, thus reducing the net hearing load and increasing the 
load anglo. The same analysis can be applied to the other 
examples although it is somewhat more complicated in the dyna- 
mic cases. Richardson did not compute the spring and damping 
coefficients and so no direct comparison of these can be made, 
but by comparing the bearing loads and load angles for the 
dynamic cases,, it can be deduced that the same order of mag- 
nitude of differences in the spring and damping coefficients 
would exist by virtue of equations k-S and 46, 

The purpose of Tables 1 and 2 is to show the small effect 
that friction has on the film forces and spring and damping co- 
efficients, and therefore to justify neglecting friction in the 
preparation of the dimensionless tables in Appendix C. In gen- 
eral, friction can be neglected in all cases except vrhen 
(CD- 20) and e are very small ^ as In half speed whirl. 

Table 3 illustrates the effect of the choice of cavitation 
boundary conditions on the spring and damping coefficients. 

Using Ockvirk’s boundary conditions. Holmes assumes that the 
entire half of the bearing is cavitated , an assijimptlon which 
approaches this theory only at the higher eccentricity ratios. 

Thus there is close agreement at n = .6, while at n = .2 the 
tvro theories yield somewhat different values. At first, it 
would be expected that this theory should alvrays give less 
values of and at low eccentricity ratios 

IuJl SS xjT xS 

since these coefficients go to zero in the completely uncavi- 
tated situation(see Appendix C), but this is not so because 
the film and the film forces obtained by this theory are un- 
symmetrical about 0 = 1.5 . and thus a clear comparison can- 
not be made. However, by comparing the results of this inves- 
tigation for a high and a lov7 feed pressure, it can be seen 
that as the feed pressure increases causing the film extent to 
increase, the radial spring coefficient decreases. 
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Sternllcht^l3j used Sommerfeld's boundary conditions and 
digital nethods to solve Reynolds’ three dimensional equation 
and to obtain spring and damping coefficients for a full Jour- 
nal bearing. However, no quantitative basis for comparison 
could be found, although it is significant to note that Stern- 
licht found the damping cross terns and 3^^. negligible, 

while this theory and Holmes’ analysis found these coefficients 
to be of the same order of magnitude as 3 and B . 

The general validity of the results can best be ascertained 
by examining the Initial assiimptions that were ma.de. Since 
low to moderately loaded bearings are considered in this investi- 
gation, bearing distortion and viscosity variations due to pres- 
sure can be neglected. Also the choice of a circumferential 
groove feeding arrangement reduces the variation of viscosity 
due to temperature since the oil flow for a line source is 
three to four times that for other feeding arrangements at any 
given load. 14 . The accuracy of the short bearing approximation 
is illustrated in Appendix A where the error is negligible for 
B/D = I', and is acceptable for L/D = |. 

The assumptions that undoubtedly need considerably more 
investigation are those that apply to the flov; of the fluid 
through the cavity, under dynamic conditions, particularly 
where the striations Join the downstream' boundary. The steady 
state investigations (referenced in the Introduction) using 
transparent bearings have shown a finite and irregular distri- 
bution of striations, with the result that considerable mixing 
must occur at the doimstream boundary, and thus in actuality 
the dovjnstream boundary cannot be precisely defined mathema- 
tically. Also, because of the loi^ pressure gradients immedi- 
ately surrounding the boundary, surface tension and vrettability 
probably influence the shape of the boimdary. 

The other assumptions that were made such as laminar flov'jr, 
Nev/tonlan fluid, etc., are generally accepted and will not be 
discussed in this paper. 
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CONCLUSIONS AND RECOMMENDATIONS 



This investigation has analyzed and programmed the un- 
steady motion of a circular, circumferentially grooved ^journal 
bearing. The results should prove useful in designing the 
particular type of bearing considered, vrhile the analytical 
and computational methods employed can be modified for appli- 
cation to other bearing configurations . The most important 
point that was brought out by this investigation was that ne- 
glecting the oil film between 6 = n and 6 = 2tt will not 
yield good values for the spring and damping coefficients in 
lightly loaded high speed bearings, where Ockvirk’s theory 
that a given load establishes the attitude angle is invalid. 
Conversely, it would be expected that assuming a complete 
film in that region would also give poor results with the 
exception of hydrostatic bearings. 

Recommendations for future work in the sub;}ect area of 
this investigation are listed below. 

1. Devise and construct an experimental apparatus to 
measure dynamic forces and load angles, vrith sufficient 
accuracy to establish the validity of the boundary 
conditions of this investigation as compared to 
Ockvlrk's boundary conditions. 

2. Extend Richardson’s theory and his program to obtain 
the spring and damping coefficients. If the differ- 
ences are relatively small, Richardson’s analysis 
would be preferable due to the simplicity of his 
boundary conditions; in its general form it could 
more readily be modified for various particular anal- 
ysis such as half speed whirl, stability, vibrations, 
etc . 

3. For use with L/D>-g, program Reynolds’ three dimensional 
equation using the boundary conditions of contiuity. 
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APPENDIX A 



Accuracy of the Short Bearing Approximation 

Donaldson[3?i obtained an exact solution to Reynolds' 
equation for a Journal bearing using series solutions and ap- 
plying Sommerfeld's boundary conditions, but he assumed that 
the feed pressure was sufficiently high to suppress cavita- 
tion. Figure XVIII is a comparison of Donaldson's solution 
and the results of this investigation, using data for uncavi- 
tated bearings from Appendix C. Two length to diameter ratios 
are plotted as a function of eccentricity ratio versus the 
ratio of approximate to exact load capacity. 

Figure XVIII 




0 0.2 0.4 0.6 0.8 1.0 



Eccentricity Ratio 
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APPENDIX B 



Computational Procedure 



The analysis was programmed in Fortran II for use on the 
IBM 709^ digital computer. The program in its present form can 
"be used either as a main program or as a subroutine by the in- 
sertion of the proper administrative and format statements. 

A sufficient number of comment statements is in tno program 
to permit someone who is familiar with the language to follow 
the logic, but some of the detailed computational methods need 
further explanation. 

It can be seen in Figure VIII that as the circumferential 
flow becomes much less than the axial flov: that the dov^nstream 
boundary conditions approach the upstream boundary conditions, 
and equation 22 applies throughout. If the flow is negative, 
l,e,, ( (A) - 20) is negative, then the positions of the upstream 
and dotijnstream boundaries are reversed. To determine which 
condition prevailed, the following expression v;as used: 



Velocity Ratio 



(VR) 



circumferential flow 


axial flow 



-A. nc 

1 Anc| +•! 6 1 



This gave not only the direction of the oirouiaferontial 
but also its relative magnitude. If VR was positive, all 
iterations were done in the positive 6 direction, and if VR 
was negative, then all iterations vxere done in the negative 
© direction. 

The beginning of the upstream boundary was obtained by 
iterating in 3 degree Intervals, starting at © = rr and calcu- 
lating w by equation 22, until the following conditions were 
agui 



1, w(I) < L 

2, w(I) < w(I -1) 
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If these conditions were not satisfied vrhen the program had 
iterated through 2tt, then a non-cavitated solution was obtained 
by means of definite integrals (see Appendix D). 

If cavitation was present, then ©2 was located by iterat- 
ing in one degree intervals starting at 0, , again calculating 
w by equation 22, until v^(l) > vr(l-l) i.e., the sign of the 
slope of the cavitation boundary changed. This established 0^ . 
One degree Intervals vxere used between 0) and ©2 so that greater 
accuracy could be obtained in calculating h* later in the pro- 
gram. 

Then, starting again at 6, or 0, +3* whichever was even, 
the program ;^as iterated in 6 degree Intervals using equation 
22 to calculate the film width along the upstream boundary, and 
equation 2? to calculate the film width along the downstream 
boundary. Equation 2? was applied as follovj-s: 

1. Assume that w(I) = w(I-l) 

2. Determine h* by iterating through the w's previously 
calculated between ©, and ©^ , until w* > w(l). 
then Interpolate between the stored values of ©* on 
either side of vj(1) to obtain h-“-. 

3. Calculate the slope ^ by means of equation 27, using 
h*(I-l) and w(I-l). 

4. Then, w(I) w(I-l) + || A© 

Hovjever, because equation 2? becomes Invalid, both mathemati- 
cally and analytically, as_A_-^0, and because the quantity 
(h-h*) becomes very small and practically impossible to calcu- 
late digitally with any accuracy as VR — ^ 0, it was necessary 
to apply certain restrictions along the downstream boundary. 
These vjere; 
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1. If 1 vr| > . 25 . equation 27 applied, 

2, If 1 vr 1<,25, equation 22 applied. 

The number .25 was obtained by trial and error; as VR app- 
roached this value, the downstream boundaries by equations 22 
and 27 agreed very closely (Figures XV and XVI have VR values 
of .33 -.33 respectively). At lower values of VR , this 

program gave inconsistant results for the reasons mentioned. 

For each © , the differential forces and coefficients 

T^ere calculated by equations 28 and 28 and equations 37 through 
44. The Integral of h*dz for use in calculating the friction- 
al forces in the cavitated zone v?as evaluated by Simpson’s Rule 
using 3 ordinates ; 

I h^dz ib ^ (h* - h*)(vr* - ^^ave 

Then the differential forces and coefficients were integrat- 
ed from © = TT to © = 3tt, by means of the subroutine .which 
uses Simpson’s Rule and 6l ordinates. For clarity, each compon- 
ent of each equation was calculated separately; thus DCSSl is 

the first term of dK„^. 

ss 

The remainder of the program is self explanatory. It 
should be noted though that FR is the radial force directed 
Inward; this was changed for conventions sake. 
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FORTRAN 

FORTRAN 

ANALYSIS CF A C IRCUNFlREN T I A^LY GROOVED JOURNAL GEARING 



INPUT IS TOT\L BEARING L ENGTh, JOURNAL RADIUS , RAU I AL CLEARAr.Cw, 
ANGULAR ^ PECi' ,GR3 IT AL SPEED, RACIAL VELOC I TY , ECCENTRI C I TY RATIO, 
ABSOLUTE FEED PRESSURE, AND VISCOSITY ( E L , R , C , OM , PH I DO T , D E , EN , PO , U ) 
UNITS ARE LO, IN, SEC, RAD IANS 

CUTPUT IS RADIAL, TANGENTIAL, HORIZONTAL , AND VERTICAL BEARING FORCES 
SPRING CCEFFF IC lENTS , and DAMPING COE F F I C I ENT S ( FR , F T . r- TOT , ANG , C RR , 
CSS,CSR,CRJ,r;KR,BSS,uRS,CXX,CYY,CYX,CXY,cXX,jYY,DXY,) 



L I MENS I Ok X2 ( 99 ) , w2 ( SD ) , X ( 6 1 ) , V; ( 6 1 ) , OF 1 1 ( 6 1 ) , DF 1 2 ( 6 1 ) , DF 1 3 ( 6 1 ) , 
1CF1A(61 ) , CF2. (6 1 ) ,OFji ( 61 ) , D F 2A ( 6 1 ) , DCRR 1 ( 61 ) ,DCRR2 ( 61 ) ,DCSRl (61 ) , 
2CCSR2 ( 61 ) ,AHSDZ ( 61 ) , DCRR3( 61 ) , DCRS4( 61 ) , DCSSA ( 61 ) , DC SR3 ( 6 1 ) , 

3CF15 (61 ) , DF2S (6 1 ) ,W1 ( 12A) 

COMBINE AND REDUCE INPUT 



GA=0M-2.0*PHID0T 
F=P0-15.0 
P I = 3 . 1 A 1 6 
EE = i.j.5*cL 
GG = ( G . D «EL ) ** 3 
CC= (U’~ R) / (2.0i^C««3) 
CD= ( U*OMsR-'^-c': )/C 
SS = Cl//C 
TT = CC->OM 



0S = C. 

E N = 0 . 

EP = 0. 

Bt = 0. 

CR = 0. 

CS = 0. 

CT = 0, 

CU^O. 

CV = 0. 

DY = C. 

AHSDZ=G. 

PAGNITUDE AND DIRECTION OF CIRCUMFERENTIAL FLOW 

VR = GA/ ( ABSF( GA) +A8SF ( DE/ ( EN^^C ) ) ) 

IF(VR) 210,110,110 



LOCATE XI FOR NEGATIVE FLOW BY ITERATING IN 3 DEG INTERVALS IN 
NEGATIVE THETA C I REG 1 I ON , B EG I NN I N G AT 5 PI 



210 CO 220 K=1,12A 
IF(K-12A) 2i:,300,999 

211 A=K 
X=K 

X 1 = B . «■ P I -i- . 05 2 35 - ( AS . 052 35 ) 
AhlS = C->( 1 .0 + EN«-CCSF ( XI ) ) 
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/iAA = COSr { XI ) 
OBB=SI \F ( XI ) 

EQUATION 20 



PH=.137t)<U=^EL*»2»(-2.CisD: »A AA + GA^» hN*-C BBci ) / ( Cc- { 1 .0 + £N»AAA ))»»•• 3 
IF{A.0*P11 + P) 212i"-:5,215 

EQUATION 22 



212 ( K) = SQ<1TF( - { ( C» ( i .0 + EN*AAA ) ) ) / { 3. C<^U* { (-2 .0--DE* 

1 AAA) + ( GAs-ENi^CaBad ) ) ) ) 

IF(K-l) 999,220,213 

213 IF(W1(K-1 )-.-99aEL) 220,214,214 

214 IF(Wi(K)-Wi{Q-l ) ) 240,220,220 

215 U1 ( K ) =0. 5*EL 
22w CONTINUE 



LOCATE X2 FOR NEGATIVE FLO'a BY ITERATING IN ONE DEG INTERVALS IN 
NEGATIVE THEIA OIRECTICN , BEGINNING AT XI 



240 CO 250 J=l,99 

241 A^J 

X2{ J)=Xl-(Aa.01745) 
rOUATIDN 22 



0^^2 { J ) = SQRTF{ - ( ( Ca { 1 .0 + ENaCCSF( X2( J )) ) ) a* 3aP ) / ( 3. 0* U* ( ( -2 . OaDE a 
1C0SF(X2( J) ) ) + {GAahNaCaSINF(X2( J) ) ) ) ) ) 

IFU-3 ) 250,250,245 
24d IF(W2{ J)-W2( J-3) ) 250,251,251 

250 CONTINUE 

251 'aQ 2 = W2 { j-3 ) 

XX2=X2 ( J-3) 

JJ=J-12 

I=61-K/2 

GO TO 153 



LOCATE Xi FO.-, POSITIVE FLOQ BY ITERATING IN 3 DEG INTERVALS IN 
POSITIVE THETA C I RE CT ! ON , 3 EG I NN ING AT 3 PI 



110 LO 120 K=l,l?4 

IF (K-124) 111,300,959 
Hi A = K 
K = K 

Xl=3. 0363+1 Aa. 05235 ) +2.0aPI 
AhlS=Ca( 1.0+ENaCaSF(Xl) ) 
AAA=COSF ( XI ) 

EQB=SINF(X1) 

EQUATION 20 



FH= . 1875 aUa£L a* 2*‘ ( ~ 2 .0*DuaAAA + GAa£NaC*8BB ) / ( C® { 1 ,0+tNaAAA ) ) a»3 
IF(4.0aph+P) 112,115,115 
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EQUATION' 22 



112 M ( K ) = SGkTF( - ( { CM 1 .C + tN^-AAA ) ) «- 3 - p ) / ( 3. C>;-U^ { (-2 .0<-DC = 
IAAA) + (GA-Ei\<fC-BB3) ) ) ) 

IP(K-l) 990, 120, 113 
i:3 IF(V^1 ( K-i )-.A99<EL) 12G,11A,114 

119 ir ( Ivl ( K ) -U1 ( K-1 ) ) 140,120, 120 

115 'rtl(K)=0.5<cL 

120 CONTINUE 



LOCATE X2 F0'< POSITIVE FLOV^ BY ITERATING IN ONE DEG INTERVALS I 
POSITIVE THETA DIRECTION , BEGINNING AT XI 



140 CO 150 J=l,99 
J = J 

141 A=J 

X2(J)=X1+(A«. 01745) 

EQUATION 22 

0^2 { J )=SQRTF( - ( { C*( 1 . C + £N*CCSF{ X2( J ) ) ) ) ** 3*P ) / ( 3 . 0-^ 0<- ( { -2 . 0*0£ * 
1C0SF(X2 { J ) ) H (GA« CN*C'-:S INF( X2( J ) ) ) ) ) ) 

IF(J-3) 150,150,145 

145 IF(W2( J)-W2(0-3) ) 150,151,151 

150 CC'MINUE 

151 VnW2 = W2{J-3) 

XX2=X2 ( J-3 ) 

JJ=J-12 

I=K/2 

IF(2*I-K) 15^,153,999 

152 1=1+1 



INCREMENT IN 6 DEGREE INTERVALS AROUND THE JOURNAL BEGINNING AT 
XKIF XI EVEN) OR AT X1 + 3IIF XI ODD) 



153 CO 170 N=l,61 
A=I 

X(I)=3.0*PI-.10 47t-(Att.i047) 

AAA = C0SF { X( I ) ) 

BB3=SINF IX( I ) ) 

iF RELATIVE MAGNITUDE OF C I RCUMFE REN f I AL FLOW IS SMALL, 
EQUATION 22 .^PPLYS ThROUGHCUT 

IF( ABSF ( VR)-.25 ) 154,159,159 

159 IF(VR) 2604,1604,1604 

2604 IF(XX2-X(D) 154,154,166 
1604 IF(XX2-X(D) 166,154,154 

EOUATION 20 

154 FH=.18 75*U*EL*«-2-{-2.0*-0£*^AAA+GA*EN«C*8BB)/( C*-(1.0 + EN«AAA) )*-3 
IF(4.0iPH+P) 165,167,167 
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LHSTR:a.v W ( OUATIC., 22) 

16p •-. ( I) = SC.<rF { - ( ( C=- ( i . G + EM*AAA ) ) ■- v 3 iP ) / ( 3 . c*U«- ( (-2 AAA ) + (oA*»EN 

I *C*0BB ) ) ) ) 

AhS = C*» ( 1 .0 + £\*AAA ) 

XS = X ( I ) 

GB TO 168 

CGW.\STR£AM V; IS COMPUTED FROM STATEMENT 166 THROUGH 1667 

166 IF(VR) 11,15,15 

11 IF (1-61) 12,13,999 

12 k(I)=W(I+l) 

GO TO 20 

13 V>(I)=W(1) 

GO TO 20 

15 IF{ I-l ) 999, 16, 17 

16 U( I )='.-.'( 61 ) 

GO TO 20 

17 U(I)=W(I-1) 

F* FOR EQUATION 27 

20 CO 1664 J=i,JJ 

IF(J-JJ) 36,165,999 
36 IF(U( I )-W2( J ) ) 1664,1661,16 61 

1561 IFU-l )999, 1662, 1663 
1662 XS=X2(J) 

GO TO 1665 

166 3 XS = X2 { J-1 ) + ( .V ( I )-w2( J-1) )*(X2(J)-X2( J-1) )/(W2(J)-W2( J-1) ) 

GO TO 1665 

1664 CONTINUE 

1665 AHS = C»(1.0 + E\'«C0SF(XS) ) 

Ah=C<-( 1.0 + ENs-AAA) 

EQUATION 27 

CQ= (0. 5*W( I ) ;:{-2.0*DE--AAA + GA*£N*C*BBB ) + ( ( { 6 . 0-U^^l. ( I ) ) ) ) / 

1 (CA^:- { AH-AHS) ) 

1667 Vn ( I ) =W ( I ) +A8SF( DW )*. 1047 



59 IF{ W( I)-EL*^0. 5) 168, 167,167 

167 K(I)=0.5kEL 

168 AA=( 1.0 + EN*A^.A) «x-3 
ec=u ( I ) ^^43 

G0= ( 1 . 0+cN*AAA> 

RR = CQ»‘:2 

Y V. = ( 1 . 0 + £ N A A A ) * - 2 

CIFFERENTIAL FUNCriONS OF FILM WIDTH AND TRIGONOMETRIC FUNCTIONS 
IN CAVITATED ZONE 

REMOVE AU. CARDS DOWN TO 1683 IF FRICTIONAL EFFECTS ARE TO BE 
NEGLECTED. 

1682 AhSDZ=0. 

IF ( .49<;EL-W( I ) ) 1683,1683,1671 

1671 AH 1 2 S = AH 1 S + ( AHS — AH 1 S ) / 2 . 0 
CDD=( { AH12S/C )-1.0)/EN 
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;<12 = AC0SF(CC.-) 

FEE-SI NF { X12 > 

E'jUATICN 20 

Hh=. 18 7 5 iU~E^-=*2« (-2.Gi:DE*CDD + GA--2\eC<EEE ) / (C^f(1.0 + 2NfiDDD) )£•!•• 3 
IF(A.C=Ph+8) 1673,1672,1672 

1672 £lE=-EEE 



EQUATION 22 



16 7 3 '/vl2 = S0RTF (-( ( C» ( 1 . O + EN^COD ) )^*3*P)/(5.0i^U*{( -2 . 0 «-D E «• ODD ) + ( GA&EN 
1 «G« cE£ ) ) ) ) 

EQUATION 22 

US = SORTF{-((C«( 1 .0 + £i\«C0SF{ XS ) ) ) * P ) / ( 3 . 0*^ U* ( ( -2 . 0*DE 
ICOSF(XS) ) + (GA«ENs-G<SINF(XS) ) ) > ) 

SSS=ABSF( (AHlS-Ah£)/2.0) 

AhSCZ = ( SSS/3 .0 )*( wl (K ) +A.0<-W12 + WS ) 

1633 CGRR3( I ) =AHSnZ*:BBB*AAA/AA 
CCRSA ( I ) =AHSDZ<^AAA/QC 
CCSS4 ( I ) -AHSOZ*BBo/CG 
CCSR3( I )=AHSnZ« AAA*»2/AA 
CF15( I )=AHSDZ«BSB/OQ 
CF25 ( I ) =AHSDZ*AAA/QO 

CIFFERENIIAL FUNGTIONS OF FILM WIDTH AND TRIGONOMETRIC FUNCTIONS 
IN COMPLETE FILM ZONE 

CFll ( I ) = ( B3«•AAA*^f2)/AA 

CF12(I )= ( BB«AAA-BBB) /AA 

E F 1 3 ( I ) = W ( I ) AA A 

CF14( n=W( I ) »OBE/{ l.C + EN*-AAA) 

CF22( I )- ( BB»B6B*«-2)/AA 

CF23( I )-W( I ) -B8B 

CF24 ( I )=W ( I ) sAAA/ ( 1 .C+EN^fAAA ) 

CCRRl ( I ) = B3«riBB*AAA»t*2/RR 
CCRR2( I ) = W( I ) ■cAAAsBab/QQ 
CCSRl ( I ) =B3«BEfc*a?«AAA/RR 
CCSR2( I ) =W( I ) i:-AAA*«2/QC 
IF( .49^>EL-W( I ) ) 1685, 16 8 5, A2C 

420 ,'' = 174+6«I 
IF(M-360) 422,422,421 

421 M=M-360 

422 ZZ = 2.0«W ( !)/i:L 

INSERT PRINT STATEMENT HERE IF THETA AND FILM WIDTH! PER CENT) IN 
CAVITATEU REGION IS CESIREC. 

1635 IF(VRi 1686,1668,1686 
1686 1=1-1 

IF(I) 999,1687,170 
1687 1=61 

GO TO 170 
1638 1=1+1 

IF(I-62) 170,1689,995 

1689 1=1 
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170 ClNTINU 



IMTEGRAT: FU.CTIQnjS of FILR V.IDTH A.\U differential TRIGONCyE fkic 
FUNCTIONS, USING SINRSCNS RULE 

tA =SI NP ( DFl : ) 

EC SINt^(QF12) 

ED •■= S INPICF 3) 

EG ^ SIMP (OF. 2) 

EH = S IMP ( OF : 3) 

OR=SIMP( CCRRl ) 
tV=SlMP( CCSRI ) 

CELETE NEXT 10 CARDS IF FRICTIONAL EFFECTS ARE TO BE NEGLECTED 

EW = SIMP(.DCSR. ) 

BS = SIMP( CCRR.. ) 
hE = SIMP(DF14) 

EP = SIMPIDF2A) 

CR=SIMP(DCRR3) 

CS=SIMP( DCRSA) 

CT=SIMP{ CCSSn ) 

CU=S!MP( DCSR3 ) 

CV=SIMP{ DF15 ) 

CY=SIMP( DF25) 

COMPONENTS OF EQUATIONS 28, 29, AND 37-44 

F 1 1 =— 2 .0 *CCii0E*BA 
F12 = CC*GA«£Ni^C*BC 
F13 = 0» 5'“'P*R*uD 
F14 = -DD^:-BE 
F15=-SS*DV 
F21=-2 .0 »CC*DEkBC 
F22 = CC*GA^^EN<C»^DG 
F23 = 0. 5*P*K^^BH 
F24=DD*BP 
F25=SS«DY 

CRRl = TTii ( 8C- i.0K£N«8R) 

CRR2=SS*BS 
CRR3=2.0«SS*I:R/C 
ERR = -2 .0«CCs-l3A 
CRS1=-TT*BG 
CRS2=-P*R-»BH/2. C«C«EN 
CRS3=-SS»8P/EN 
CKS4 = -SS«DS/C^>EN 
ERS = -2 .0->CC«-3C 
CSS1=TT«:.C 

CSS2 = P<R«BC/? -O^C^^EN 

CSS3=-SS -BE/.N 

CSS4=-SS*0r/C*^EN 

E SS = -2 . 0 *CCK-hG 

CSR1 = TT* (BG-3.0<fEN*BV) 

CSR2=-SS»BW 

CSR3=-2.0<fSS<DU/C 

ESR=-2 .o«cc*ac 



GO TO 350 
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FORCES AKD Si^RING A\C CAMPI,\G CONSTANTS IF NC CAVITATION 



300 PRINT 912 

AI=-PI »LN/( 1 .0-lN<s 2 )<s2.5 

AG= ( 3. 1A/6N) • (1 .0-{ 1 .0/ ( 1 .C-EN<-«2 ) <« 0.5 ) ) 

AC = i.57/{1.0-EN<-c2)« = 1.5 

F 1 1=0. 

F12=-0.0 
F13 = -0 .0 
F 1A=0. 

F 15=0. 

F21=0. 

F22=CC*GA<GG“EN«2 .0*AC*C 

F23=0.0 

F24=0. 

F25=0. 

CK1U=0. 

CKR2=0. 

CRR3=0.0 

EKR=0.0 

CRS1=-GG*TT*2.0«AC 

CRS2=0.0 

CRS3=-(SS*EE/EN )*2.0»AG 

CRS4 = 0 .0 

eRS=0.0 

CSS1=0.0 

CSS2=0.0 

CSS3 = 0 .0 

CSSA=0.0 

ESS = -GG»4.0<-CC*AC 
CSR1=GG«-TT*( AC-3.0i:-ENsAI ) 

CSR2=0.0 

CSR3=0.0 

ESR=0.0 



350 FR=-2.0* ( F11+F12+F13+F14+F15 ) /(7. 5<EL*R) 

FT = 2 .0*!( F21 + F22+F23 + F24 + F25)/( 7.5«EL*R) 
FTOT=SQRTF(F i *«2 + FR=f*2) 

FHI=ATANF(FT/FR ) 

ANG=57 .3«PHI 

CRR=-2 .0 2-(CRR1 + CRR2 + CRR3 ) «C / ( 7 . 5 EL * R ) 
CRS=-2 .0* (CRS1+CRS2+CRS3+CRS4 ) *C / ( 7 . 5« E L*R ) 
CSS=-2 .0>(CSS1+CSS2 + CSS3 + CSSA ) «C / ( 7 . 5 « E L *R ) 
CSR = 2. 0* ( CSR.+CSR2 + CSR3 ) «-C/ ( 7. 5*EU^R ) 

E RR = -2 . 0 * BRR C OM C / ( 7 . 5 EL * R ) 

0SS=-2 .OcCMcC/( 7. 5cELcR ) 

BRS=2. OcBRScOMcC/ (7. icELcR) 

AAA=COSF ( PHI ) c*2 
E133 = SI NF ( PHI ) cc2 
CCC=C0SF(PHI )cSINF(PhI) 

CXX=CRR< AAA+ ( CSR-CRS ) cCCC+CSScCBB 
CYY = CRRc-BBB+ ( CRS-CSR } cCCC + CSScAAA 
CXY = CRS*AAA-i- ( CRR-CSS ) cCCC + CSRcBBB 
CYX = CRSc5B3+( CSS-CRR ) «CCC-s-CSRcAAA 
EXX = 3RRcAAA+ ■'.Oc&RS-cCCC + BSScBBB 
E Y Y = B R R c 8 B B - ; . 0 -c B R S c C C C + 3 S S c A A A 
EXY = BRRcCCC-ERScAAA-8SSa-CCC+BRScB88 



999 CALL EXIT 
END 
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F CRT RAN 



SIMPSONS RUL FOR 61 CRCINATES, KITH 6 DEGREE INTERVALS 

FUNCTION SIMMY) 

CIMEMSIGn Y(61) 

CJD=0. 

EV£N = 0 . 

DO 800 1^2,60,2 
800 EV£N = EVEi\ + Y( i ) 

DU 310 1=3,59,2 
810 CDD=ODC+Y(I) 

SIMP=(.1047/3.0)«(Y(i)+4.0«EVEN + 2.0<f0DD + Y{61)) 

RET U RN 

END 

DATA 
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APPEMDIX C 



Tabulated Values of Bearing Load, Attitude Angle, 

And Spring and Damping Constants 

Tables ^ through 45 contain values of dimensionless load, 
load angle, and dimensionless spring and damping constants as a 
function of eccentricity ratio, pressure ratio, and dimension- 
less squeeze. The eccentricity is in Increments of 0.1 from 
0.1 to 0.7, the latter being chosen as the upper limit in 
order to remain v:ithln the limits of the short bearing approx- 
imation, the assumption of constant viscosity, and the assump- 
tion of no journal or bearing distortion. The pressure ratios 
are 1.5, 2.0, 3»0. 4.0, 5»0, and 6.0, vrhich adequately cover 
the range of most bearing applications. The choice for varia- 
tion of squeeze was made as follows: from Shavr and Hacks 

p q 

small bearings : ~ & 

c .003 

P "1 

large bearings: ^ ^ 

Now let ^ Vary from ^ to 1. Then R varies from 2L to §L 

and; 






r 



^L! 



hU)i 

2 



12 ( .003)^ 

r 2 > 



I 2 

= .278 X lO" 



! c i .0005 



= 1.6 X lo’' 



Now assume a viscosity range from 10~“ to 10 Ib-sec/in^ and 
a speed range from 100 to 2000 rad/sec. Then: 






/^— A ^ /L \ 3 

Pa 
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.185 




( 10 "^’ ) ( 100 ) ( . 278 ) ( 10 ^ ) 
15 




(10"^ ) (2000) (1.6) do'* ) 
15 



21400 



However, no bearing applications at the upper limit could bo 
conceived and so it was decided to vary VV'-' from .1 to 8000 
in 30 steps that would permit easy interpolation. It should 
be noted that generally the first few entries in each table 
are uncavitated bearings, which can be recognized by virtue of 
the fact that 0 = 90^* . 

The tabulations are in terms of the entire bearing length 
(2L) and are converted to dimensional values as follows: 



F = F p^LR 



K 






KNH 



c 






(nrm 



BNM 



61? c 

p^lr 



For end fed bearings of total length L, divide the above 
values by 2. 
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TABLE h ni VENSICNLESS BEARING LOADS, LOAD ANGLES, AND SPRING AND DAMPING CONSTANTS 
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APPENDIX D 



Table of Integrals 



The following definite Integrals are used In evaluating 
equations 28 and 29 and equations 37 through 44 for the non- 
cavltated case. 




o 




0 




27T 

sin^ e de 



(1 +n cos 0) 



o 




2TT 

cos* 0 d0 



(1 n cos 6)* 



0 




b 



6 



I 



2jr 

sin 6 cos* 0 do 
(1 + n cos 0)'^ 



0 
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